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Utilizing the generalized factorization of twist-4 processes, we derive the modified heavy quark frag- 
mentation function after considering the gluon radiation induced by multiple scattering in DIS. It 
is found that the mass effects of heavy quark may reduce the gluon formation time and change the 
medium size dependence of heavy quark energy loss. The radiative energy loss is also significantly 
suppressed relative to a light quark due to the dead-cone effect. 



In relativistic heavy-ion collisions, when 
a fast parton propagating in a dense medium, 
it may experience multiple scattering with 
other partons in nucleus and lose a large 
amount of energy loss via induced gluon 
radiation 1 . Such kind of jet quenching phe- 
nomenon has aroused a lot of attention 
recently 2 ' 3,4,5 , and it is found that the total 
energy loss of a massless parton (light quark 
or gluon) has a quadratic dependence on the 
medium size due to non-Abelian Landau- 
Pomeranchuk-Migdal (LPM) interference ef- 
fect. In order to get a complete understand- 
ing of the mechanism of jet quenching it is im- 
portant to investigate the heavy quark energy 
loss by multiple scattering in nuclei with uni- 
form methods. Currently, though the frame- 
work of the light quark energy loss induced 
by gluon radiation has been well established, 
there are only few studies about the energy 
loss of a heavy quark in literature 6 > 7 ' 8 - 9 . Here 
I would like to report our group's study on 
jet quenching of massive quark in nuclei with 
the twist expansion approach 8 by utilizing 
the generalized factorization in pQCD 8,10 . 

To separate the complication of heavy 
quark production and propagation, we con- 
sider a simple process of charm quark pro- 
duction via the charge-current interaction in 
DIS off a large nucleus. The results can be 
easily extended to heavy quark propagation 
in other dense media. 

To the leading-twist in collinear approx- 
imation, the semi-inclusive cross section fac- 



torizes into the product of quark distribution 
fs e ( x b+xm), the heavy quark fragmentation 
function Dq^h(zh) (%h — ^h/^q) an d the 
hard partonic part HjS)(k, q, M) 8 
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As illustrated in Fig. 1, in medium, the 
propagating heavy quark in DIS will suffer 
additional scattering with other partons from 
the nucleus, which will induce gluon radia- 
tion and cause the leading quark to lose en- 
ergy. Such induced gluon radiations will ef- 
fectively give rise to additional terms in the 
evolution equation leading to the modifica- 
tion of the heavy quark fragmentation func- 
tions in nuclei. There are total 23 cut dia- 
grams that contribute to these double scat- 
tering processes. 

Similar to the case of light quark prop- 
agation in nuclear medium 10 , the gener- 
alized factorization of twist-4 processes will 
be employed to evaluate these contributions. 
The dominant contribution gives the semi- 
inclusive tensor for heavy quark fragmenta- 
tion from double quark-gluon scattering, 
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Figure 1. A typical diagram for double scattering 
processes with three possible cuts. 



+ (g — frag.) + (virtual corrections) . 

The above is very similar to the case of dou- 
ble scattering of a light quark 10 and resem- 
bles that of gluon radiation off a heavy quark 
in vacuum. The transverse momentum dis- 
tribution + (1 - z) 2 M 2 ] 2 is typical 
of bremsstrahlung from a heavy particle. It 
vanishes in the small angle It — > 0. One 
can rewrite such a suppression relative to the 
gluon radiation from a light quark as 
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where 9o = M/q~~ and 9 = lx/q~z. This 
is often referred to as the "dead-cone" phe- 
nomenon that suppresses small angle gluon 
radiation and therefore reduces radiative en- 
ergy loss of a heavy quark 6 . 

The contribution of double scattering is 
proportional to a twist-4 parton correlation 
function 
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in the limit xl *C x 8 ' 10 , where xl = 
£^/2p + q~ z(\ — z), xa = 1/ttinRa- The 
coefficients a\ and a 2 are polynomial func- 
tions of M 2 jl\ and become (1 + z)/2 and 
C F (1 - z) 2 /2C A , respectively for M = 8 . 

Rewriting the sum of single and double 
scattering contributions in a factorized form 
for the semi-inclusive hadronic tensor, one 
can define a modified effective fragmentation 
function Dq^h(zh, At 2 ) as 

Dq^h(zh,V 2 ) = D Q ^ H {z H ,y 2 ) 
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where D Q ^ >H (z H ,y 2 ) and D g ^ H (z H , fi 2 ) are 
the leading-twist fragmentation functions of 
the heavy quark and gluon. The modified 
splitting functions are given as 
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and A^ q ^ gq (z) = A^ q ^ qg (l - z). 

The gluon formation time for radiation 
from a heavy quark can be read out from the 
phase factors in the effective twist-four ma- 
trix element as 

1 2z(l 
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p+x L £ 2 T + {1- z) 2 M 2 ' 

which is shorter than that for gluon radiation 
from a light quark. This should have signif- 
icant consequences for the effective modified 
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quark fragmentation function and the heavy 
quark energy loss. 

One can then calculate the heavy quark 
energy loss, defined as the fractional energy 
carried by the radiated gluon, 
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where xm = (1 — z)M 2 /2zp + q~ and x^ = 
/j, 2 /2p + q~ z(l — z) + xm- Note that xl/xa = 
L a /t/ with L~ A = Ration /p + the nuclear size 
in the chosen frame. The LPM interference 
is clearly contained in the second term of the 
integrand that has a suppression factor 1 — 

~2 / 2 
£ X L/ X A, 

bmce xl jx A ~ x B M 2 /x A Q 2 , there are 
two distinct limiting behaviors of the energy 
loss for different values of xb/Q 2 relative 
to x A /M 2 . When x B /Q 2 > x A /M 2 for 
small quark energy (large Xb) or small Q 2 , 
the formation time of gluon radiation off a 
heavy quark is always smaller than the nu- 
clear size. In this case, 1 — exp (—x L /x A ) ~ 1, 
so that there is no destructive LPM interfer- 
ence. The integral in Eq. (6) is independent 
of Ra, and the heavy quark energy loss 
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is linear in nuclear size Ra- In the oppo- 
site limit, xb/Q 2 *C xa/M 2 , for large quark 
energy (small xb) or large Q 2 , the quark 
mass becomes negligible. The gluon forma- 
tion time could still be much larger than 
the nuclear size. The LPM suppression fac- 
tor 1 — cxp (— x\ /x 2 A ) will limit the available 
phase space for gluon radiation. Therefore, 
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Figure 2. The nuclear size, Ra, dependence of charm 
quark energy loss for different values of Q 2 and xg. 
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Figure 3. The Q 2 dependence of the ratio between 
charm quark and light quark energy loss in a large 
nucleus. 



the heavy quark energy loss 
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now has a quadratic dependence on the nu- 
clear size similar to the light quark en- 
ergy loss. Shown in Fig. 2 arc the 
numerical results of the Ra dependence 
of charm quark energy loss, rescaled by 
C(Q 2 )C A a 2 (Q 2 )/N c , for different values of 
xb and Q 2 . One can clearly see that the Ra 
dependence is quadratic for large values of 
Q 2 or small xb- The dependence becomes 
almost linear for small Q 2 or large xb- Here 
we take M = 1.5 GeV for charm quark in the 
numerical calculation. 

To illustrate the mass suppression of ra- 
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diative energy loss imposed by the "dead- 

cone , we plot the ratio -^r ot 

charm quark and light quark energy loss as 
functions of Q 2 and Bjorkcn variable in 
Fig. 3 and Fig. 4 respectively. Apparently, 
the heavy quark energy loss induced by gluon 
radiation is significantly suppressed as com- 
pared to a light quark when xb is large and 
when the momentum scale Q, or the quark 
initial energy q~ is not too large as com- 
pared to the quark mass. Only in the limit 
M <C Q, q~, is the mass effect negligible. 
Then the energy loss approaches that of a 
light quark. 

In summary, we have calculated the en- 
ergy loss of the massive quark in terms of 
modified heavy quark fragmentation function 
in nuclei with the twist expansion approach. 
We show that mass effects such as 'dead- 
cone" effect gives a significant suppression to 
the induced heavy quark energy loss. In par- 
ticular, (nuclear) medium size dependence of 
heavy quark energy loss is found to change 
from a linear to a quadratic form when the 
initial energy and momentum scale arc in- 
creasing, which is quite different from the 
light quark energy loss where the total en- 
ergy loss always have a quadratic dependence 
on medium size. The result can be easily ex- 
tended to a hot and dense medium, which 
may be applied to study heavy quark produc- 
tion and suppression in heavy-ion collisions at 
RHIC and LHC. 
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Figure 4. The xb dependence of the ratio between 
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